Osmotic compression of droplets of hard rods: 
A computer simulation study 



Yu. Trukhina,^ S. Jungblut,^ P. van dcr Schoot,'^ and T. Schilling'^ 

^ Institut fiir Physik, Johannes Gutenberg- Universitat, D-55099 Mainz, Staudinger Weg 7, Germany 
^Institut fur Physik, Universitat Wien, Boltzmanngasse 5, 1090 Wien, Austria 
'^Faculteit Techmsche Natuurkunde, Technische Universiteit Eindhoven, 
Postbus 513, 5600 MB Eindhoven, The Netherlands 
(Dated: February 12, 2009) 

By means of computer simulations we study how droplets of hard, rod- like particles optimize 
their shape and internal structure under the influence of the osmotic compression caused by the 
presence of spherical particles that act as depletion agents. At sufficiently high osmotic pressures 
the rods that make up the drops spontaneously align to turn them into uniaxial nematic liquid 
crystalline droplets. The nematic droplets or "tactoids" that are formed this way are not spherical 
but elongated, resulting from the competition between the anisotropic surface tension and the 
elastic deformation of the director field. In agreement with recent theoretical predictions we find 
that sufficiently small tactoids have a uniform director field, whilst large ones are characterized by a 
bipolar director field. From the shape and director-field transformation of the droplets we are able 
to estimate the surface anchoring strength and an average of the elastic constants of the hard-rod 
nematic. 

PACS numbers: 61.20. Ja,64.70.Md,64.70.Ja 

I. INTRODUCTION 

Fluids of elongated particles that interact via a harshly repulsive excluded- volume interaction potential, i.e., hard- 
rod fluids, have proven immensely useful as a model system for investigating the formation, structure and properties 
of liquid crystals!. Indeed, hard-rod fluids exhibit a wealth of aggregated states including nematic, smectic and 
columnar liquid- crystalline phases as well as various plastic and crystalline solid phases, depending on the precise 
shape, density and composition of the model particles present in the system?. Closest experimental realizations of 
the hard-rod fluid model are found in fluid dispersions of very stiff polymers, inorganic rod-like colloids, filamentous 
viruses, fibrillar or tubular protein assemblies, and carbon nanotubes^ii^i^i^ii. For this class of system, the transition 
from an isotropic to a symmetry-broken liquid-crystalline or crystalline state is entropy rather than enthalpy driven, 
and the relevant control parameter is not the temperature but the concentration or particle density. The transitions 
between the different phases are driven by a tradeoff between free volume and configurational entropy. 

Arguably, the experimentally and theoretically by far most extensively studied liquid-crystalline state is the nematic 
that is characterized by uniaxial and up-down symmetry. The preferred mean orientation of the particles is described 
by a director field the ground state of which is spatially uniform. The transition from the isotropic to the nematic 
phase is first order. (See, however, the work by Oakes and co-workers^.) In the biphasic region it proceeds via 
the formation of nematic droplets called tactoids in the background isotropic dispersion. Tactoids are typically not 
spherical but elongated and spindle-shaped, and they have over the past century or so been observed in a wide variety 
of systema^ i-'^^i-'^-'^i-'^^'-'^'^i-'^^i-'^^i-'^^ . The shape and internal structure of nematic droplets is the result of the competition 
between the preferred surface anchoring of the director field and the deformation of the director field that occurs if 
the preferred anchoring is indeed accommodated. 

The shape and director-field configurations of tactoids have recently been studied theoretically in considerable 
detail within a macroscopic Frank elasticity theorjsi^iiiii^iiSii^. Predictions depend on several parameters: two surface 
energies (the surface tension and anchoring energy), three bulk elastic constants associated with the splay, twist and 
bend deformations, the saddle-splay surface elastic constant, and, finally, the size of the droplet. For our purposes, 
it suffices to summarize the most important predictions, presuming preferential planar anchoring of the director field 
to the interface that for rod-like particles seems to hold for entropy reasons: i) The director field of small drops is 
uniform and of large ones bipolar; ii) The crossover is smooth and set by a healing or extrapolation length that is the 
ratio of an average elastic constant and an anchoring energy; iii) The aspect ratio of uniform tactoids depends only 
on the anchoring energy, and of bipolar ones on their size relative to the healing length. 

Some aspects of these theoretical predictions have been verified against experiment, in particular the size dependence 
of aspect ratio and opening angle of the sharp ends of the tactoid a^'^i^^'^^i"'^^'^" . The crossover from bipolar to uniform 
director field has not been observed and this probably presents quite an experimental challenge because it has been 
predicted to occur when the drops are quite small, that is, in the micrometer range. For such small droplets, not at 
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all very much larger than the particles that they are made up from, director field patterns are difficult to distinguish 
in polarization microscopic images. 

Clearly, computer simulations are of use her o^^'^^i^^ , not least because a macroscopic description might break down 
for such small drops, in other words, the predicted transition may be spurious. On the other hand, recent simulations 
on the nucleation of the nematic phase in a hard-rod fluid have indeed shown that small nuclei of the nematic phase 
are elongated and it appears at least that their director field is uniform^^. Simulations on small nematic droplets 
of prolate particles interacting via a Gay-Berne potential have been shown to exhibit similar behavior, although 
fluctuations are large so a director field is not so easily definedi^i^ Because systematic computational investigations 
of the shape and internal structure of nematic droplets are lacking, we set out to do this for nematics of hard 
spherocylinders of aspect ratio eleven. Our simulations build on our earlier work on hard-rod nematics in a spherical 
cavity of fixed size and shape^l, but we now focus on actual drops that can adjust their structure. This we do by 
placing the spherocylinders in a suspension of spherical particles, which interact via a hard-core repulsion with the 
spherocylinders yet are transparent with respect to themselves, i.e., behave like an ideal gas. 

The method allows us to compress by osmotic stress droplets of a fixed and predestined number of rod-like particles. 
This is achieved in a controlled manner because the density of the rods in the drop will be set by the imposed pressure 
of the ideal gas of spheres (if we ignore a small contribution from a Laplace pressure). Hence, the drop size must then 
be given by the number of rods in the simulation box. Note that our making use of what in essence is a depletion 
agent (the spherical particles) is not all that dissimilar in spirit to the experiments of Dogic and Fraden^-, who used 
the non- adsorbing polymer dextran to induce phase separation in aqueous dispersions of filamentous fd virus particles. 

Our results may be summarized as follows. 

i) For relatively low densities of spheres, the droplets are spherical and the rods randomly oriented; 

ii) Above a critical sphere density, the drops become nematic and elongated; 

iii) Small factoids have a uniform director field and large ones a bipolar one; 

iv) At the crossover the drops measure about four rod lengths so indeed are very small; 

v) The surface tension anisotropy and healing length that we deduce from the simulations are consistent with 
theoretical estimates. 

The remainder of this paper is organized as follows. In section [TT] we introduce the model and define the relevant 
order parameters. A detailed discussion of our results is given in section [TTTl Finally, in section HVl we summarize our 
findings and relate them to theory and experiment. 

II. MODEL AND SIMULATION DETAILS 

In the Monte-Carlo simulations, we model the rods as spherocylinders consisting of a cylindrical part of length 
L and diameter D, capped at both ends by hemispheres of diameter D. In order to stabilize droplets of rods, we 
surround them by a liquid of spheres. The spherocylinders interact via a hard-core repulsion. The spheres also interact 
with a hard-core repulsion with spherocylinders but are able to freely interpenetrate each other at not energy cost. 
Hence, our model is an extension of the Asakura-Oosawa or AO model for colloid-polymer mixture a^^'^^ to anisotropic 
colloids^°. Note that strictly speaking the spheres do interact with each other via the rods whilst this is not so in 
the original AO model. This difference is of academic interest only, because we consider only those concentrations of 
particles that produce a very strong phase separation into sphere- and rod-rich phases, so the former behaves as a 
(nearly) ideal gas of spheres that in essence acts as a barostat for the droplet consisting of only spherocylinders. 

As already advertised in the introduction, this implies that by means of changing the number of spheres in the 
system at fixed volume, we can tune the pressure exerted on the rods in the droplet. Through that we vary the 
density of the drop and therefore also the physical properties of the drop, such as the elastic response if the drop is 
in a nematic fluid state. Hence, we use the number density of spheres psp measured far away from the droplet as a 
parameter that characterizes the external conditions imposed on the droplets formed. The same method has been 
used in an earlier study to study the formation of nanocrystals in a simulations^. 

The simulations were performed at fixed particle number N and simulation box volume V (and temperature T, but 
as the system is purely entropic, temperature is not relevant here) in a cubic box with periodic boundary conditions. 
We focus on spherocylinders with an aspect ratio of_L/£'-t-l = Il (i.e. L/D = 10), and spheres with a diameter 
twice that of the thickness of the rods, i^sp = 2_D. Spheres that are much larger than this introduce strong effects on 
the surface anchoring, whilst for smaller spheres the numbers needed in the simulation are so large to be unpractical 
from a computational point of view. 
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The number of spherocylinders in the box was varied from 200 to 700, and the number of spheres was fixed such that 
phase separation was induced into two phases containing virtually only spherocylinders or spheres. This corresponded 
in our simulations to between approximately 20000 and 70000 spheres. We chose the simulation box to be 70^D'^, 
i.e., sufficiently large to ensure that the spherocylinders did not interact directly with each other via the periodic 
boundaries. We verified the droplets that form are not system spanning. The systems were equilibrated by local 
translation and rotation moves. Depending on the specific concentrations, 10^ — 10^ MC sweeps were required. 

The boundary and volume of a droplet is established as follows. We first divide the system into small boxes and 
next verify whether a box contains the center of a sphere or whether it is intersected by any spherocylinder. In the 
latter case, this box is counted as a part of the droplet. If it so happens that a box does not fall into either of the two 
categories, it is counted as a part of the droplet if it has a larger number of nearest neighbouring boxes containing 
spherocylinders than spheres. 

As is customary, we define the average alignment of the (N) rods in terms of the traceless tensor Q, with the 
elements 

1 ^ 

Qal3 = ^ ^ {iUiaUip - Sap) , 
i=l 

where Uia is the a component (a = x,y,z) of the unit vector along the axis of particle i = 1,...,N and 5a[3 the 
Kronecker delta. Diagonalization of the tensor yields three eigenvalues, A+, Aq and A_, where A+ > Aq > A_. Different 
authors use different combinations of these eigenvalues to define the nematic and biaxial order parameter o^^'^'^ . To 
avoid confusion, we shall present all eigenvalues instead. The case A+ > 0, Aq = A_ corresponds to a structure with 
one preferred direction. The case A+ = Aq > corresponds to a structure in which one direction is avoided and the 
two other directions are equally favored. All the intermediate cases A+ > 0, A-|_ > Ao > A_ correspond to a biaxial 
structure. Obviously, in an isotropic phase one finds A+ = Aq = A_ = 0. 

When calculating observables, such as the components of the order parameter tensor properly averaged over the 
ensemble of configurations, all the possible types of symmetry in the system have to be taken into account. In the 
case of radial symmetry we calculated the orientational tensor directly by averaging over all configurations obtained 
in the simulations. 

If there is a preferred axis in the system along which the particles tend to align, we need to proceed differently. 
As rotations of the director do not cost any (free) energy, the orientation of this axis can fluctuate strongly during 
a simulation run. In order to average the local properties of interest, the configurations have been rotated in such a 
way that the director always points in the same direction. It is important to point out that this procedure introduces 
additional noise to the simulation results. In the following, we show all data rotated such that the director is aligned 
with the z-axis. 



III. RESULTS AND DISCUSSION 

In the diagram of fig. [1] we have indicated with asterisks the system conditions for which we ran the simulations. 
The letters "a" through "1" are used in the following to refer to specific points in this diagram. We have indicated 
schematically in the figure the shape and the structure of the droplets, which we discuss in more detail below. The 
numbers in the boxes indicate the aspect ratios of the droplets. 

Figs. [2] and [3] give the snapshots of two typical structures given by conditions "i" and "c", respectively, where only 
the spherocylinders are shown for clarity. Two essential differences between these structures can be seen: The first of 
the two droplets shown is i) more elongated and ii) the spherocylinders in this droplet have a much stronger tendency 
to orient in one direction. Below we will analyse these effects in a more quantitative fashion. 

The difference is caused by the difference in bulk sphere density and hence osmotic stress imposed on the drop by the 
sphere fluid. Droplet "i" is subject to much higher pressure and hence is condensed much more than droplet "c" , and 
has crossed over to the nematic phase. The rod densities pD^ in the two drops are 0.034 and 0.027 (averaged over each 
droplet). From the (remarkably accurate) Lee-Parsons theory of the nematic transition in bulk fluids of hard rods^, 
we expect for spherocylinders of aspect ratio 11 the nematic transition to occur at a pressure of PD^/ksT « 0.25, 
corresponding to an ideal gas density of PspD^ ~ 0.25. This is in reasonable agreement with the transitional regime 
around PspD^ w 0.23 in the diagram. We also estimated the transition density by simulating a compression curve 
and an expansion curve in the bulk and obtained PspD^ « 0.22, which again is in good agreement with the numbers 
above. We note that due to the effects of Laplace pressure the nematic transition should occur at a somewhat higher 
sphere density then that in bulk solution. For details we refer to the Appendix. 

In flg. m we have plotted a typical iso-density distribution in the r-z-plane for nematic droplet "i" . Inside the 
droplet the spherocylinders have an approximately constant rod density equal to pD^ = 0.039, which rapidly decays 
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when approaching the interface to the fluid of spheres. It is clear that the drop is not spherical, the aspect ratio 
being approximately 1.8. For comparison, and in order to determine the shape of the droplets, we cut a slice from 
the r — 2-density profile at half of its maximum value. That is where we expect the Gibbs plane to be situated. The 
curves obtained in this way are shown in fig. [5] for several systems consisting of the same amount or spherocylinders 
(^rods = 500) but with various densities of spheres, corresponding to the points "a" , "b" , "c" , "f " , "i" , and "k". 

At low pressures the droplets are (nearly) spherical. This is to be expected if the rods are in their isotropic state. 
However, if the pressure is increased the droplets crossover to the nematic phase and hence become elongated in order 
to reduce either the elastic deformation of the director field if the anchoring is strong, or the anchoring free energy 
if the anchoring is weak. As already advertised, the crossover occurs when the typical drop dimension exceeds the 
healing length, which we are going to discuss in more depth in the following section. 

On increase of the density of rods, the elastic constants in all likelihood increase, too, as should the interfacial 
tension and potentially also the anchoring strength. We expect from scaling arguments that the ratio of the anchoring 
strength and the surface tension is a weak function of the pressure^, however, and in the weak anchoring regime it is 
this ratio that dictates the aspect ratio of the drops^^. In the strong anchoring regime the aspect ratio is an increasing 
function of the healing length that in the equal-constant approximation is given by the ratio of the elastic constant 
and the surface tension. Hence, whether the droplets become more elongated then depends on how strongly these two 
energies depend on the pressure. Apparently, the elastic constants increase more strongly with pressure because the 
aspect ratio increases from about 1.1 to 1.8 with increasing sphere concentration. 

We observe the same tendency for the systems of 700 spherocylinders. At low pressures the droplet is more or less 
spherical (state point "e", aspect ratio 1.1) and at high pressures the droplet becomes elongated (state point "j", 
aspect ratio 1.8). We have to note, however, that droplet size also affects the aspect ratio: it decreases with increasing 
size. This has been observed experimentally and is predicted theoretically based on macroscopic theoryi^ii^ii^. The 
aspect ratios for all the state points are shown in the boxes in the "shape" diagram fig. [1] with an estimated error of 
0.1. 

We now turn to the orientational state of the rods in the droplets. This is described in terms of the eigenvalues 
of the orientational tensor Q (see|n|. In fig. [5] we show the r- and z-profiles of these eigenvalues for the drops "i", 
"f", and "c" (that each consist of 500 spherocylinders). The droplets become less ordered and less elongated. Droplet 
"i" has a high nematic ordering in its center, which slightly decays on approach of the interface to the gas of spheres 
in both, the r- and the z-direction. Droplet "f" has a lower density and also a nematic order parameter that is 
smaller. Droplet "c" has even lower values of the nematic order parameter, and at the interfaces all the eigenvalues 
become equal to 0, indicating an isotropic configuration. This drop probably is close to the conditions where the 
isotropic-nematic transition takes place. 

By analysing the tensor Q for all the investigated droplets, we construct scalar order parameter profiles as well as 
the nematic director-field configuration. The results are indicated schematically in the diagram [TJ 

1) Droplets "a" and "b" are spherical droplets of an isotropic rod fluid; 

2) Droplets "c", "d" and "e" are in the isotropic-nematic transition region, exhibiting a strongly fluctuating 
orientational order; 

3) States "f" to "1" are strongly nematic drops; 

4) The director field of drops "i", "j", "k" and "1" is more or less uniform, those of "g" and"h" bipolar, and of "f" 
in between these two. 

To illustrate these findings, we show the director field (the axis given by the eigenvector corresponding to the largest 
eigenvalue of the orientation tensor) as a function of the radial and axial distances r and z in fig. [7| for the case "k" 
and in fig. [5] for the case "h" . The director field of the factoid "k" is more or less uniform and oriented along the 
main axis of the drop. That of "h" is bipolar, i.e., curved along the elongated drop surface toward the tips, where the 
scalar order parameter drops to zero. This signifies the melting of the nematic near the tips, where theoretically one 
would expect the surface point defects ("boojums") to reside^^. 

The director field in a nematic droplet is determined by the interplay of surface anchoring and elastic forces. The 
bipolar structure can only be formed if the energy for bending is small enough compared to the surface energy. This 
can be achieved if the density of the suspension is small (but the droplet still has a nematic structure), or if the 
droplet is big (and, therefore, the curvature of the interface is small). Systems "g" and "h" fit into this category. 

The aspect ratio of the droplets does not depend on their size if the director field inside the droplets is homogeneous. 
This follows directly from the well-known Wulff construction of the droplet shape given any polar angle-dependent 
surface tensionii. Indeed, all the droplets with an approximately homogeneous director field, droplets "i" , "j" , "k" , 
and "1" , that we obtained in our simulations have the same aspect ratio of 1.8 within a statistical error of the simulation 
of about 0.1, in agreement with this theoretical prediction. From the aspect ratio observed in the simulations we 
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can in fact deduce a dimensionless anchoring strength. Let us presume that the anisotropic surface tension 7 has 
a functional form of the Rapini-Papoular type, so 7 = t[1 + ijj(q • n)^], with r the bare surface tension, uj the 
dimensionless anchoring strength, q the surface normal and n the director field at the surface of the drop. For planar 
anchoring to be favored, uj > 0. From the Wulff construction wc then deduce that the aspect ratio of the drop equals 
1 + (jj for < w < 1 and for oj > 1^. Hence, we find from our simulations a value for uj of 0.8, quite close the 
value of 0.65 found by Dijkstra and co-workers in a simulation study of nematic drops nucleated in a super-saturated 
dispersion of hard rods^^. 

Finally, from the crossover from bipolar to homogeneous director fields, we obtain an estimate for the healing length. 
According to fig. 1, the crossover occurs for drops of a volume about 1.7 x 10'*£''^. This value may actually depend on 
the sphere density, but because of the lack of any detailed information we shall ignore this for simplicity. Presuming 
that the bend elastic constant is about ten times larger than the splay elastic constant, which seems reasonable on 
account of predictions for hard rods in the Onsager limit, the crossover occurs at a droplet volume equal to about 
10 times the healing length A = {Ku — K24)/tu! cubed according to macroscopic theor y^^i^^ , where Ku denotes the 
splay elastic constant and K24 the saddle-splay surface elastic constant. So, we find for the healing length A w 12D, 
which is about a rod length. Clearly, macroscopic theories, such as those of Kaznacheev and collaboratorsl^°, and 
of Prinsen and van der Schoo t^^i^^'^^ , could perhaps be expected to break down at such small length scales, yet the 
predicted crossover from uniform to bipolar director fields apparently still survives. 

From the estimate of the healing length, we can obtain an order of magnitude estimate of the interfacial tension 
between the rods and the spheres that we can compare with the scaling estimate given in the Appendix. If we presume 
the Saupe-Nehring relation to hold, implying that K24 — {Ku — K22)/2 — , and make use of the approximate 
expression K22 = ifii/3 we obtain A « 2^11/3™. Hence, Ku/t « 141?, or, i3tD^ « f3KiiD/U. Within 
a second-virial approximation, which admittedly is not very accurate for rods of aspect ratio below 20, we expect 
f3KiiD Ri 0.9 to hold near the transition^^. Hence, for the interfacial tension between the rods and the spheres we 
obtain the estimate firD^ k, 0.07. According to our scaling estimate cited in the Appendix, we have firD^ k, 0.25a^/Z? 
with a a prefactor that should be of order 0.1— and ^ the interfacial width. If ^ « L, this then implies that a « 0.03, 
which is somewhat smaller than expected^. 

In order to go beyond this quantitative analysis, simulation data of the elastic constants of the bulk nematic and of 
the surface tension between the co-existing bulk fiuids would be necessary. Unfortunately, these are not available yet 
for spherocylinders of aspect ratio 11, as, in particular, simulations to determine elastic constants are computationally 
rather expensive. (To our knowledge, elastic constants have been computed only for spherocylinders of aspect ratio 



By means of computer simulations we have shown that fluid droplets of hard rods, osmotically compressed by the 
presence of spherical particles, undergo an isotropic-nematic transition at sufficiently high osmotic stress. We find the 
nematic droplets not to be spherical but elongated. The director field of the drops is uniform if smaller than a critical 
size and bipolar if larger than that. We interpret our findings in terms of the predictions of continuum mechanical 
theory that minimizes the combined effect of an elastic deformation of the director field and an anchoring frustration 
of this director field at the surface of the drops. Although in our simulations the drops are not at all large on the 
scale of the rods, and continuum theory should perhaps not be expected to be accurate, results from both levels of 
description seem to be consistent with each other down to drop sizes that are as small as a few times the particle 
length. 



An estimate for the rod density in the droplet as a function of the number N^-ods of rods and the density of spheres 
Psp can be obtained by presuming complete demixing of the two components and by presuming that the interface 
between them is sharp. Let R be the radius of the drop, assumed perfectly spherical, then p — 3A'rods/47ri?^ is the 
density of the rods in drop. The bulk pressure of the hard rods is to a very good approximation equal to the expression 
put forward by Parsons and by Lee^ 



IV. SUMMARY 



V. 



APPENDIX 




at least in the isotropic phase, where (j) = p{7tD^ /6+ttLD^ /A) denotes the packing fraction, A = {L / D)"^ /tt{1 + 3L /2D) 
for slender rods is proportional to their aspect ratio and /3 denotes the reciprocal thermal energy l/ksT with ks 
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Boltzmann's constant and T the absolute temperature. The pressure of the ideal gas of spheres obeys 

Mechanical equilibrium between the gas of hard rods and that of ideal spheres demands that 

27 

p J L — p 

R ~ 

where the second term on the left-hand side is the contribution from the Laplace pressure across the curved interface, 
with 7 the interfacial tension that presumably depends on the bulk densities of both the rods and the spheres. For 
any given number of rods N^^dsi this equation sets the equilibrium size of the drop. 

An estimate of the magnitude of the Laplace pressure may be given by making use of the scaling Ansatz 7 « Psp£,j 
with ^ w L the actual interfacial widths. Hence, P/Psp « 1 — 2a£,/R with a a constant of proportionality that we 
estimate to be of order 0.1 . So, the presence of the interface reduces the pressure in the drop relative to that in 
the reservoir of spheres and hence postpones the onset of the nematic phase to higher densities of spheres the smaller 
the drop. 

From Table I of Lee^, we deduce by linear extrapolation that for rods of L/D ~ 10 the bulk nematic phase sets 
in at a dimcnsionless pressure /3PD^ « 0.247, corresponding to a sphere fraction of PspD^ w 0.247, in reasonable 
agreement with what we find in the simulations. See fig. 1. 
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FIG. 1: "Shape" diagram of hard sphcrocylindor droplets immersed in a fluid of spherical particles. The horizontal axis shows 
the volume of the droplet of spherocylindcrs in units of cylinder thickness cubed, D'^ , and the vertical axis the number density of 
the spheres far away from the droplet. The letters "a" to "1" are used in the main text to refer to specific points on the diagram. 
The sketches show schematically the shape and the interal structure of the droplets, and distinguish spherical isotropic droplets, 
elongated nematic droplets with either a homogeneous or a bipolar director field. The numbers in the boxes indicate the aspect 
ratios of the droplets. The boundaries of the boxes distinguish between: isotropic (dashed), transition region (dash-dotted), 
bipolar (dotted) and homogeneous (solid). 
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FIG. 2: A snapshot of a nematic droplet of rods that forms for the conditions indicated by "i" in fig. 1. The spherical particles 
are not shown for clarity. The nematic total director is along the main axis of the drop. The average rod density in the drop 
is equal to pD^ = 0.034. 
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FIG. 3: A snapshot of an almost isotropic droplet of rods for conditions "c". The average rod density in the drop is equal to 
pD^ = 0.027. 
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FIG. 4: Isodensity lines of spherocylinders in a droplet from conditions "i", shown in cylindrical coordinates. The solid line 
marks the boundary of the drop. Successive dashed, dash-dotted, etc., line demarcate densities from pD^ = 0.039 to pD^ = 0. 
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FIG. 5: The shapes of droplets consisting of 500 sphcrocylinders at different pressures of spheres (conditions "a", "b", "c", "f", 
"i", and "k" indicated in the schematic of fig. 1). Droplets "a" and "b" arc those of isotropically oriented rods, and hence 



more or less spherical, 
nematic drops, with "i 



Drop "c" is in the transition zone from the isotropic to the nematic state. Drops "f" 
' and "k" exhibiting more or less uniform director fields and "f" a more bipolar one. . 



"i" and "k" are 



13 



r [D] 



z[D] 




system c 

( ^o o o o o o o o o e-^.^^ 



15 



10 

z[D] 



1 

0.75 
0.5 
0.25 


-0.25 



FIG. 6: The r- and z-profiles of the eigenvalues of the orientational tensor Q for the systems "i", "f, and "c" (each consisting 
of 500 spherocyUnders). 
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FIG. 7: An example of a nematic droplet with a homogeneous director field (system "k"), the eigenvectors corresponging to 
the maximum eigenvalues of the orientational tensor Q are shown in polar coordinates. 
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FIG. 8: An example of a nematic droplet with a bipolar director field (system "h"), the eigenvectors corresponging to 
maximum eigenvalues of the orientational tensor Q are shown in polar coordinates. 



